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Abstract. The mean area of a two-dimensional Gaussian ring of N monomers is known to diverge when 
the ring is subject to a critical pressure differential, p c ~ N^ 1 . In a recent publication [Eur. Phys. J. E 
19, 461 (2006)] we have shown that for an inextensible freely jointed ring this divergence turns into a 
second-order transition from a crumpled state, where the mean area scales as (A) ~ N, to a smooth state 
with (A) ~ N 2 . In the current work we extend these two models to the case where the swelling of the ring 
is caused by trapped ideal-gas particles. The Gaussian model is solved exactly, and the freely jointed one 
is treated using a Flory argument, mean-field theory, and Monte Carlo simulations. For a fixed number Q 
of trapped particles the criticality disappears in both models through an unusual mechanism, arising from 
the absence of an area constraint. In the Gaussian case the ring swells to such a mean area, {A) ~ NQ, 
that the pressure exerted by the particles is at p c for any Q. In the freely jointed model the mean area 
is such that the particle pressure is always higher than p c , and (A) consequently follows a single scaling 
law, (A) ~ N 2 f(Q/N), for any Q. By contrast, when the particles are in contact with a reservoir of fixed 
chemical potential, the criticality is retained. Thus, the two ensembles are manifestly inequivalent in these 
systems. 

PACS. 36.20.Ey Macromolecules and polymer molecules: Conformation (statistics and dynamics) - 
64.60.Cn Order-disorder transformations; statistical mechanics of model systems - 61.25.Hq Swelling 
of polymers 



1 Introduction 

Polymer rings confined to two dimensions (2D) have been 
the subject of several theoretical works |TJ 2,3,4.5, 6]ITI[S1 
[ ^ flCl l [TT 1 [T2 l fl^ l fT4]. in part as highly idealized models for 
membrane vesicles. The statistics of 2D chains have been 
studied also experimentally, using polymers adsorbed at 
liquid interfaces [15|I16| or membranes [17] . as well as vi- 
brated granular chains [TB] and rings [TH] • In those experi- 
mental systems both random-walk and self-avoiding-walk 
statistics were observed. 

Pressurized 2D rings were theoretically investigated in 
a number of works p[7l[8PIT0lfllTIT2lfl^^ . In the Rudnick- 
Gaspari model PfTOlfTT] the ring is represented by a closed 
Gaussian chain of N springs subject to an inflating pres- 
sure differential p. This model yields a divergent mean area 
at a critical pressure, p c ~ iV -1 . Recently we have demon- 
strated that the inclusion of chain inextensibility (i.e., 
considering a freely jointed chain rather than a Gaussian 
one) turns this divergence into a second-order smoothen- 
ing transition [14j . For p < p c the ring is in a crum- 
pled, random-walk state with the mean area scaling as 
(A) ~ Nf < (pN), whereas for p > p c the ring is smooth 
with (A) ~ N 2 f > (pN). 

In many practical cases a pressure difference across a 
closed surface is achieved osmotically through a concen- 
tration difference, e.g., by trapping particles inside the 



envelope. In the current study we extend the models of 
pressurized Gaussian and freely jointed 2D rings to the 
case where the swelling is caused by a trapped ideal gas 
of particles. We examine both the canonical and grand- 
canonical ensembles, in which the particle number Q or 
the particle fugacity 7 are fixed, respectively. 

The model is defined in the following section. In Sec- 
tion [3] we derive exact results for a Gaussian ring swollen 
by trapped particles. We then treat the swelling of a freely 
jointed ring in Section [4] by employing a Flory argument, 
mean-field calculation and Monte Carlo simulations. In 
addition, we derive exact asymptotes for the large-swelling 
regime. Finally, the results are discussed in Section [5] 



2 Model 

We model a 2D polymer ring as a closed ideal chain of 
N monomers and no bending rigidity. The ring encapsu- 
lates an ideal gas of particles. We shall consider separately 
the two cases of a Gaussian and a freely jointed chain. 
The unperturbed monomer-monomer link length is taken 
as the unit length, 1 = 1. For a freely jointed chain this 
length remains fixed under swelling. For a Gaussian chain, 
however, the link length has a statistical distribution that 
varies with swelling; / is then defined as the unperturbed 
root-mean-square link length. 
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In a canonical ensemble, where the number Q of par- 
ticles is fixed, the partition function is given, up to a con- 
stant prefactor, by 



Z(N,Q)= dAP (N,A)A Q /Q\, 



(1) 



where A is the area bounded by the ring, and Pq(N,A) 
is the probability distribution function of the area for an 
unpressurized TV-monomer ring (for either a Gaussian or 
a freely jointed chain). The mean area of the ring is given 

by 

(A(N, Q)> = (Q + 1)Z(N, Q + l)/Z(N, Q). (2) 

In a grand-canonical ensemble, where the particles are 
in contact with a reservoir of fixed fugacity 7, we use equa- 
tion ([T]) to write the grand partition function as 



Z(N,j) 



00 />oo 

V Z(N, Q) 7 « = / dAP (N, A)e~< A . 

Q=0 J ° 



(3) 



The mean area in this ensemble is given, therefore, by 

(A(jV, 7 )) =ainz/a 7 , (4) 

and the mean particle number by 

<Q(iV,7)) = dlnZ/dln-y = 1 {A). (5) 

Thus, the mean particle density, c = (Q) / (A), is equal to 
the fugacity 7, as it should for an ideal gas. 

The correlation between fluctuations in particle num- 
ber and in the ring area is characterized by the following 
covariance, 



C QA {N n ) 



(QA)-(Q)(A) 

(QKA) ■ 



From equation ([3]) we find 



(QA) 



1 d 2 z 

Z (9(ln7)c>7 



= l(A 2 ), 



which, combined with equation (0, yields 



Cqa(N, 7 ) 



(AA 2 ) 



(6) 



(7) 



(8) 



where (A A 2 ) = (A 2 ) — (A) 2 is the mean-square area fluc- 
tuation. Hence, interestingly, the cross-correlation of Q 
and A is identical to the relative mean-square fluctuation 
in the ring area. 

Our model system consists of two subsystems — parti- 
cles and ring monomers — which may have different relax- 
ation times. We note, however, that in equations ([1]) and 

all configurations of both subsystems are traced over. 
Hence, the equilibrium properties studied in this paper are 
unaffected by such time-scale differences. 

Several complications related to the definition of the 
ring area in the model should be mentioned. In principle 



the A appearing in equation |T]) should be the geomet- 
rical area A g of the ring, since it is A g that determines 
the translational entropy of the particles. As in previous 
works .!) 10 miTTI . however, we are technically bound to 
use the algebraic area A instead, i.e., the area obtained 
from an integral over the ring contour. This area may 
contain both positive and negative contributions, A + and 
In an unperturbed ring positive and negative areas 
are equally favorable, and the mean algebraic area then 
vanishes. (The algebraic area may also "count" a certain 
geometrical-area contribution more than once due to chain 
winding.) On the one hand, particle entropy forces us to 
use only the positive-area part of Pq (hence the integra- 
tion from to 00 in equation ([1])). Consequently, the mean 
algebraic area of an unperturbed ring, (A), does not van- 
ish. This is the definition used in the analytical parts of 
this work. In the simulations, on the other hand, we follow 
changes in the actual algebraic area of the ring, allowing 
the total area to become negative. Thus, in the numerical 
parts of this work (A) of an unperturbed ring does vanish. 
In addition, the simulated particles are placed only inside 
positive parts of the ring area, i.e., their entropy-relevant 
area is A + ^ A. Finally, in equation ([3]) one notices a di- 
rect analogy between the problem of fixed particle fugacity 
studied here and that of an empty ring subject to a fixed 
pressure, which was studied in [9jfT4] . The mapping is not 
exact, however, since in the current case A is restricted to 
positive values, whereas in [9lH"4] it was not. All these sub- 
tleties are significant only in the weak-swelling regime of 
small areas. As the inner pressure exerted by the particles 
increases, the distinction between the various area defi- 
nitions becomes negligible and does not affect our main 
results as will be demonstrated below, as will be demon- 
strated below. 



3 Gaussian Ring 

In this section we consider a 2D Gaussian ring swollen by 
trapped particles, for which exact results can be derived. 
The chain consists of a set of N springs with fixed spring 
constant, A = 1, in units of k^T /I 2 (to yield an unper- 
turbed root-mean-square spring length of I = 1). For the 
sake of the following sections, in which the spring constant 
is allowed to change, we keep the results dependent on A 
without substituting A = 1. 



3.1 Canonical Ensemble 

The probability distribution function of the algebraic area 
for a 2D Gaussian ring was calculated by Khandekar and 
Wiegel 0[3]0 



P G (iV,A) = 



1 



2ttA 



X N Ncosh 2 (2ir\A/N)' 



(9) 



1 Equation @ differs from the original Khandekar- Wiegel 
expression by a factor of 2, since the particle translational en- 
tropy in the current model restricts our analysis to positive 
areas (Eq. JTJ|). 
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Substituting equation (|9|) in equation |T]), we obtain the 
exact partition function, 

2 / N \ Q 

Z G (N,Q) = ^ J (1-2 1 -0)C(Q), (10) 

where £ is the Riemann zeta function. The mean area is 
then given, according to equation ([2]), by 



(A G (N,Q)) 



N(Q + 1) 1- 2-Q C(Q + i) 



4ttA 1 - 2 l -Q ((Q) 
which turns, in the limit of large Q, to 

(A G (N,Q > 1))=NQ/(4ttX). 



(11) 



(12) 



In the current case of fixed A, we thus get that the ring 
swells gradually (linearly) with Q. This is qualitatively 
different from the Rudnick-Gaspari result for the swelling 
as a function of fixed pressure [9l fT0lfTT] . where the mean 
area diverges at a finite critical pressure, p c = Att/N (in 
units of U-qT /I 2 ). The ability of the Gaussian ring to swell 
indefinitely stems, as in the Rudnick-Gaspari model, from 
its extensibility. 

Since the mean area is proportional to Q, we obtain in 
this ensemble the peculiar result that the particle density 
is independent of particle number, 



c G (N) =Q/(A G ) = 4nX/N. 



(13) 



The same holds for the particle fugacity, 

7 G (7V) = exp(-<9 In Z G /dQ) = AttX/N, (14) 



which is equal to the density of equation (|13|) . as expected 
for an ideal gas. The mean pressure of the gas, in our 
dimensionless units, is p = c = 4tt/N. (Recall that A = 1.) 
Thus, for any Q the ring swells to such an extent that the 
pressure exerted on it by the gas is always at the Rudnick- 
Gaspari p c , yet without the associated criticality. 



3.2 Grand-Canonical Ensemble 

The grand partition function is obtained from equations 
© and © as 

7) - [~ A + M- x ) ~ ~ x)] , (15) 

where x = jN/ (87rA), and ip is the digamma function (the 
logarithmic derivative of the gamma function). The mean 
area is calculated according to equation ([4]) , yielding 

(A Q {W)=^±M-*) (16) 
-xty'(-x)-ip'(l/2-x)]}, 

where prime denotes a first derivative. The mean number 
of trapped particles, according to equation ([S]), is simply 
given by (Q G ) = j{A G ). 



The digamma function and its derivative diverge for 
small arguments as ip(x) ~ —x~ 1 and ip'( x ) — x ~ 2 \ re ~ 
spectively. Thus, the expressions for Z G , (A G ), and {Q G ) 
diverge at a critical fugacity (with A = 1), 



7c = 4tt/7V. 



(17) 



Since for the ideal gas p = c = 7, this divergence is analo- 
gous to that of the Rudnick-Gaspari model at p c = 4n/N. 



4 Freely Jointed Ring 

We now turn to the case of a freely jointed chain, i.e., an 
inextensible ring whose link lengths are fixed. 



4.1 Canonical Ensemble 

4.1.1 Flory Argument 

We begin by examining a simple Flory argument, following 
the lines of [14]. The free energy of the ring (in units of 
kftT) is expressed as a function of i?, the radius of the 
statistical cloud of monomers (i.e., the mean radius of 
gyration). We divide it into three terms, 



F(R, Q) — F e \ + -Fincxt + Fen, 

F el ~ R 2 /N, F inext - RyN 3 , 



(18) 

Q[ln(Q/R 2 ) - 1]. 



The first term is the usual entropic-spring free energy of 
a Gaussian chain [5D]. The second is the leading non- 
Gaussian correction due to inextensibility (see App. A of 
|14j). The last term comes from the translational entropy 
of the ideal gas, where the mean area of the ring is taken 
proportional to R 2 [U[5] . 

Unlike the case of fixed pressure [14] , the free energy 
of equation (|18j) exhibits no phase transition with increas- 
ing Q. Upon minimization with respect to R we get the 
following scaling law: 



R 2 - N 2 f(Q/N), 



(19) 



which, despite the crudeness of the Flory argument, turns 
out to be correct, as demonstrated below. 



4.1.2 Mean-Field Approximation 

We proceed to study the freely jointed model by relaxing 
the strict constraints of fixed link lengths into harmonic 
potentials. This kind of approximation, used in several 
previous polymer theories 21 . 22 . 23lll4| . is equivalent to a 
mean- field assumption [22] . The calculation is identical to 
that for a Gaussian chain (Sec. [3]), yet the spring constant 
A is not fixed but determined self-consistently to maintain 
the root-mean-square link length equal to / = 1. 

To impose the relaxed inextensibility constraint we dif- 
ferentiate Z G with respect to A to get the mean-square link 
length, and then set it to 1, 



N- 1 d\nZ G /d\ = (1 + Q)/A = 1, 



(20) 
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where Q = Q /N is the rescaled particle number. For Q — > 
we get A = 1, as expected, but as the number of trapped 
particles increases, the springs get stiffer according to 

A=l + Q. (21) 

This value is substituted in equation (|12p to obtain 

(A FJ (N, Q)) = A max Q/(Q + 1). (22) 

As Q — y oo the mean area appropriately tends to its max- 
imum value, A max = N 2 /(Ait), which is the area of a circle 
with perimeter N. Equation f2"2"]) shows that the swelling 
of the ring with increasing Q is gradual, without any phase 
transition. 

The particle density is given by 



4.1.4 Monte Carlo Simulations 

We performed Monte Carlo simulations to obtain the mean 
area (A) of a freely jointed ring as a function of the num- 
ber Q of trapped ideal-gas particles for various ring sizes 
N. The algorithm combines an off-lattice scheme for the 
polymer, as was used in [14], with a lattice model for the 
particles. The simulated system is schematically shown in 
Figure [TJ The ring is represented by a polygon of N equal 
edges of length 1 = 1. The 2D coordinates of the vertices 
take continuous values. The coordinates of the particles 
are defined on a 2D square lattice with lattice constant of 
either I or 1/5, depending on the required precision. 



c F \N, Q) = Q/ (A F,] ) = 4ir(Q + 1)/N. (23) 
The fugacity is found from equations (fT4l) and (|2"T|) as 

7 FJ (^V,g) = 47r(Q + l)/iV, (24) 

which is equal to c FJ , as it should for an ideal gas. The 
exerted pressure is p = 7 = (4tt/N)(Q+1), which is larger 
than p c — 4tt/N for any Q. 



4.1.3 Exact Asymptote for Strong Swelling 

In the limit of large particle number, Q > 1, the parti- 
tion function can be calculated exactly. In this limit the 
chain is highly swollen and its statistics is governed by 
large-area conformations. Thus, to apply equation |T]), we 
merely need to know how the area probability distribution 
function for a 2D freely jointed ring, P FJ (N, A), decays to 
zero as A approaches A max . This calculation is presented 
in Appendix A. The result is 

P FJ (N, A < A niax ) ~ (A max - A) N ' 2 . (25) 

Substituting equation (|25p in equation |T]), we get 

where r is the gamma function. The mean area is then 
obtained from equation ([2]) as 

(27) 

Thus, the exact approach of the mean area to its maximum 
value differs from the mean-field result, equation (f2"2"]) . by 
a factor of 1/2. 




Fig. 1. Schematic view of a section of the simulated system. 
Point particles (marked by +) are located on lattice sites con- 
fined inside the ring. On each step each particle is moved to 
a randomly chosen neighboring site. In the second part of the 
step a single randomly chosen chain vertex (marked by a cir- 
cle) is moved to the only other position (marked by a diamond) 
that maintains the lengths of the two links attached to it. 

The initial configuration is a fully stretched, regular 
polygon of purely positive algebraic area, A = A + . Dur- 
ing the simulation we keep track of positive and negative 
contributions to A and place particles only on lattice sites 
belonging to A + . Each step of the simulation is composed 
of two actions. In the first, each particle is moved to a 
randomly chosen juxtaposed lattice site, unless this site 
is outside the positive-area part of the polygon. In the 
second action, a randomly chosen chain vertex is moved 
to the only other position that satisfies the edge-length 
constraint (see Fig. Q}. This action is automatically re- 
jected if it makes a particle leave the positive-area part 
of the ring. These dynamics involve 0(Q) operations per 
step. The number of operations required for equilibration 
is 0(N 3 Q), limiting our investigation to N < 500. The 
simulations were performed for N between 50 and 400 
and Q between and AN. 

Figure [2] shows the simulation results for the mean area 
(scaled by A max ~ N 2 ) as a function of the rescaled par- 
ticle number Q = Q/N. All data for different values of N 
collapse onto a single universal curve, in accord with the 
scaling law obtained from the Flory argument and mean- 
field calculation (Eqs. ([19]) and (|22]l). Yet, only for small Q 
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does this curve coincide with the mean-field scaling func- 
tion. For such small Q it coincides also with the Gaussian 
result, (Eq. (12]) with A = 1). 
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Fig. 2. (Color online) Mean area of a freely jointed ring, scaled 
by A max ~ iV 2 , as a function of the rescaled number of parti- 
cles. Monte Carlo simulation results for various values of N are 
represented by symbols with error bars. The dashed line shows 
the mean- field result (Eq. (J22J ) . The result for a Gaussian ring 
(Eq. (|12|) ') is given for reference (dotted line). The inset focuses 
on the small Q region using a logarithmic scale. 



To confirm the highly swollen behavior derived in Sec- 
tion 14.1.31 we simulated rings of N — 50 with Q between 
and 50AT. The results are presented in Figure El show- 
ing good agreement for large Q with the exact asymptote, 
equation (|27|) . 



4.2 Grand-Canonical Ensemble 



4.2.1 Flory Argument 



To account for contact of the particles with a reservoir 
of fixed fugacity 7 we add a — Qln 7 term to the Flory 
free energy, equation (|18|). Minimization with respect to 
Q yields the following grand potential: 



F a ~ R 2 /N, F inoxt ~ R 4 /N\ 



(28) 



- 7 i? 2 . 



This Landau-like free energy has a second-order phase 
transition at 7c ~ Af _1 , similar to the one obtained for 
fixed pressure in |14) . For 7 < 7 C Flncxt is negligible, and 
R has a Gaussian distribution with (R 2 ) ~ N(l — 7/ 7c ) _1 - 
For 7 > 7 C we have R 2 ~ N 2 (j/j c — 1). 



4.2.2 Mean-Field Approximation 

As in section Fi.l.2l we apply a relaxed inextensibility con- 
straint by determining the spring constant A so that the 
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Fig. 3. Mean area of a freely jointed ring, scaled by A max , as a 
function of the rescaled particle number, as obtained by Monte 
Carlo simulations for N = 50 (symbols with error bars). The 
solid line shows the exact asymptote for a highly swollen ring 
(Eq. (J27J ) . The dashed line presents the mean-field result (Eq. 

d). 



root-mean-square link length should be I — 1. In the cur- 
rent grand-canonical case this is done by differentiating 
Z of equation (JT5|) with respect to A, 



-N- 1 dlnZ G /dX 

1 2 
= A + \ N + 1 NZ G 



x{tp{-x) - ip(l/2 - x) 



(29) 



thus obtaining a transcendental equation for A (7, N). Equa- 
tion (|29[) can be combined with equation (fll)]) to get a 
simpler expression for (A) as a function of A, 



(^ J (iV, 7 ,A))=7V(A-l)/ 7 . 



(30) 



Numerical solution of equation (|29|) for A and substitution 
of the result in equation (|3H)l yield the mean area as a 
function of 7 and N. 

We can get a good approximation for A(7, N) by sub- 
stituting for the diverging terms in equation ([29]) ip(x) ~ 



-x 1 and ip'{x) ~ x 2 . This gives 



7 + 1 + 



jj (7 + 1) + -w 



A( 7 ,AT»1) 

(31) 

where 7 = 7 A r /(47r) is the rescaled fugacity. This result 
for A is the same as the one obtained for fixed pressure p 
14J with p = 7. In the thermodynamic limit, defined here 
as N — » 00 and 7^0 such that 7 is finite, equation (f3Tj) 
reduces to the continuous but nonanalytic function, 



A( 7 , — > 00) 




(32) 
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Substituting equation (|3"T|) in equation ([30| yields an 
approximate expression for {A} as a function of 7 and N, 

<A FJ ( 7 <1,7V»1))^ 

^ 7-l + ^ + y/(7-l) 2 + ^(7 + l) + i^ ( 33 ) 
which in the thermodynamic limit becomes 



{47r 7(1—7) 4-rr 1— 7 ' ^ 

^ | 7 -l|<JV-Va 
£^^£(7-1) 7-l»iV- 1/2 , 

(34) 

thus exhibiting a continuous phase transition at % = 1, 
analogous to the one found for fixed pressure (Eq. (16) 
in [14]). As shown in equation l|5|), the mean number of 
trapped particles is given by (Q FJ ) = r y(A FJ ) and, there- 
fore, undergoes the same transition at %. 

4.2.3 Exact Asymptote for Strong Swelling 

As demonstrated above, and also confirmed by simula- 
tions in the next section, fixing the fugacity in the grand- 
canonical ensemble is analogous to fixing the pressure in 
an empty ring. Thus, to get the asymptotic swelling at 
high fugacity we can readily use the exact asymptote for 
high pressure derived in [H] and replace p — ► 7. This 
yields 



(A( 7 >l,iV))~A : 




where I n are modified Bessel functions of the first kind. 



4.2.4 Monte Carlo Simulations 

We performed grand-canonical Monte Carlo simulations 
to check the results of the preceding sections. To each step 
of the algorithm described in section 14.1.41 a third action 
is added as follows. Either a randomly chosen particle is 
removed with probability Q/(^fA + ), or a particle is added 
to a randomly chosen site within the positive-area part of 
the ring with probability ~fA + /{Q + 1). 

Figured] shows simulation results for different values of 
N in the crumpled and smooth states, 7 < 1 and 7 > 1, re- 
spectively. The data collapse onto a single universal curve 
in each state. In the weak-swelling state (Fig. |4j\) the 
mean area scales as (A) = Nf < (^/), whereas in the strong- 
swelling state (Fig. |3p) it scales as (A) = N 2 ^^). These 
results confirm the scaling laws predicted by the Flory ar- 
gument and mean-field theory, albeit with different scal- 
ing functions. The discrepancy in Figure 2]A between the 
simulation results and the mean-field expression stems 
from the different definitions of the ring area discussed 
in section In the simulation we measure the actual al- 
gebraic area, which vanishes at 7 — 0. The calculation, 



however, considers only positive algebraic areas (cf. Eq. 
([T])) and, hence, yields a finite mean area at zero fugac- 
ity, (A(0)) = AT(ln2)/(27r). When this value of (A(Q)) is 
subtracted from the mean-field result (dotted line in Fig. 
2]^), the agreement with the simulations is excellent, indi- 
cating that the mean-field theory does accurately capture 
the swelling for 7 < 7 C - On the other hand, the disagree- 
ment between simulation and theory in Figure 2J3 arises 
from the inadequacy of the mean-field approximation for 
stretched freely jointed chains, as was observed also in 
P3]. In FigureQp we compare the data for a freely jointed 
chain subject to a pressure p > p c P3] with those for fixed 
fugacity 7 > 7 C - The two data sets are practically indis- 
tinguishable once one identifies p with 7. 

Due to computation limitations we have not directly 
confirmed the asymptotic strong-swelling behavior at large 
7 as given by equation (|35[) ; the equivalent asymptote for 
high pressure was confirmed in [14] . 

To test equation © for the cross-correlation of parti- 
cle number and ring area we measured from the simula- 
tions the covariance C QA = ((Q A) - (Q) (A)) / ((Q) (A)) for 
N = 50 and varying 7. The results are presented in Fig- 
ure [5] alongside the relative mean-square area fluctuations. 
The two data sets are indistinguishable, in agreement with 
equation ([8]). (In these measurements the area was taken 
as A+.) As expected, the correlation is appreciable in the 
crumpled, unpressurized state (small 7) and decays to zero 
as the ring swells into a smooth circle (large 7). 

5 Discussion 

Fluctuating empty envelopes, such as the polymer rings 
studied previously ^Vgy71I5UgyTDyTTlfT^llT5VTt] , have the spe- 
cial property that their thermodynamic state can be de- 
fined by specifying the surface degrees of freedom and in- 
tensive parameters only. For example, given the number of 
ring monomers N, temperature T, and pressure p of such 
an empty ring (actually, N and the ratio p/T suffice), one 
can completely characterize the system, including its ex- 
tensive mean area {A} . This situation is somewhat similar 
to the one encountered in the thermodynamics of electro- 
dynamic modes in a cavity. Since the number of photons is 
unconstrained, the system can be defined without specify- 
ing an extensive parameter; given the surface properties of 
the cavity and the radiation pressure, one can determine 
the cavity volume. 

In the current work the rings are not empty but en- 
close a gas of particles. Yet, we do not constrain their area, 
neither directly nor via a conjugate pressure. Hence, these 
systems still belong to the special category just described. 
When the number Q of trapped particles is prescribed, 
we find that the critical behaviors, previously reported 
for pressurized rings, disappear. The mechanism is qual- 
itatively different from ordinary ones, where criticality is 
removed by disorder or fluctuations. It lies here in the ab- 
sence of an area constraint, i.e., in the freedom of these 
systems to select their mean area so as to maximize the 
total entropy. For a Gaussian ring we have found that the 
mean area is (A) = NQ / (47r) , implying that the particles 
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Fig. 4. (Color online) Mean area of a freely jointed ring as 
a function of particle fugacity below (A) and above (B) the 
critical point. The fugacity is scaled by -y c = 4n/N, and the 
area by N in (A) and by v4 max = N 2 /(4n) in (B). Symbols with 
error bars show the results of grand-canonical Monte Carlo 
simulations for different values of ring sizes N. The dashed 
lines show the mean- field prediction in the limit N — > oo. The 
dotted line in (A) presents the mean-field result shifted down 
by (A(0))/N = (ln2)/(2vr). Panel (B) shows also simulation 
results from [14] for the swelling of a freely jointed ring with 
pressure p, for which the horizontal axis represents p = p/p c . 



exert an effective pressure just equal to p c — A-k/N for 
any Q > 0. The area never diverges, as it does for fixed p, 
but grows gradually, indefinitely, with increasing Q. For a 
freely jointed ring, by contrast, the mean area is such that 
the exerted pressure remains above p c for any Q > 0. As 
a result, the swollen ring always obeys the smooth-state 
scaling law, (A) ~ N 2 f(Q/N). 

When the particle fugacity 7 is prescribed, the critical- 
ity of the two models is retained. In this case the system 
does not have the freedom to select the pressure exerted on 
its boundary, since p is determined by 7 through the equa- 
tion of state of the gas (p = 7 for an ideal gas). The mod- 
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Fig. 5. Correlation between particle number and ring area as 
a function of particle fugacity, compared to the relative mean- 
square area fluctuations. Simulations were performed for an 
iV = 50 ring, and the ring area was taken as A+. 



els then become analogous to the fixed-pressure ones, as 
studied in [9|I14] . Thus, the canonical and grand-canonical 
ensembles are inequivalent in these systems. A vivid ex- 
ample is seen in the Gaussian model, where the canonical 
mean fugacity 7 is found to be independent of Q (Eq. 
(|14p ). whereas the grand-canonical mean particle number 
does vary with 7 (Eqs. |5]| and (dH)) and even diverges at 
7 C . The former result of a fugacity independent of parti- 
cle number is clearly a consequence of the unconstrained 
area. We infer, therefore, that this unusual property also 
underlies the ensemble inequivalence. 

In the grand-canonical case we have derived, and con- 
firmed in simulations, an identity relating the correlation 
between particle number and ring area with the area fluc- 
tuations (Eq. © and Fig. [5]). The only assumption enter- 
ing this derivation is the consideration of the particles as 
an ideal gas. Thus, the same identity should hold for any 
fluctuating envelope with unconstrained volume, enclosing 
an ideal gas of particles. 

This study has been restricted to self- intersecting rings. 
When a self-avoidance term is added to the Flory grand 
potential of section l4.2.1l the second-order transition dis- 
appears, as in the case of fixed pressure [13]. We note, 
however, that the Flory argument does not reproduce the 
correct scaling regimes for pressurized self-avoiding rings 
as derived in [SIEKS]. 

The pressurized 2D envelopes studied here and in [14] 
are evidently idealized systems. Possible implications for 
more realistic systems of unconstrained volume, e.g., 3D 
self- avoiding vesicles |24j with pores or aquaporin chan- 
nels, are currently under study. 

This work was supported by the US-Israel Binational Sci- 
ence Foundation (Grant no. 2002-271). 
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Appendix A 

Area probability distribution function of a highly 
swollen freely jointed ring 

Unlike the probability distribution function of the area for a 
Gaussian ring, P<p(iV, A) of equation ([9]), its counterpart for a 
freely jointed ring, Pq J (N,A), is not known analytically. The 
asymptote of Pq as A — > i majI , nonetheless, can be calculated, 
which is the aim of this Appendix. 

A configuration of the ring is defined by a set of N 2D 
vectors specifying the monomer positions, {rj}j = o....,iv with 
i"o = rjv to make the chain closed. The partition function for 
a fixed area A is given by 

N 

Z F \N,A) = J J J <lr„y. r v - l)S(A'[{rj}] - A) 



dp 



3 = 1 



ip(A'[{ rj }]-A) 



Ti-l|-l), 

(A.l) 



where ^'[{r,-}] is the area of the configuration {r?}. 

When the ring statistics is governed by highly swollen con- 
figurations, the integration over {r,} can be performed ana- 
lytically using the transfer-matrix technique |14| . This leads 
to 



Z FJ (iV,^<^ max ) 



dpe 



— ipA 



27T/ 



(ipN\ 



(A.2) 



where 7o is the zeroth-order modified Bessel function of the first 
kind. The asymptotic expansion of To for large arguments, 



1/2 



/ r .Y 

e i" . 



(A.3) 



Substituting it in equation (|A.2|) . we get, up to a constant 
prefactor, 



Z F3 (N,A< A ma x) 



/ 



dpe 



ip(A u 



-A) 



(pN) 



-N/2 



(A.4) 



which, upon a simple change of variables (or, alternatively, a 
stationary-phase approximation) readily yields 



Z tJ (N,A< A max )~(A n 



A) 



N/2 



(A.5) 



Thus, the area probability distribution function, which is 
proportional to Z FJ , is given to leading order in (j4 max — A) 
by ti 

P FJ (iV, A < A max ) ~ (A max - A) N/2 . (A.6) 

This result is used in section 14.1.3 1 to calculate the mean area 
of a freely jointed ring in the strong-swelling regime. 
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